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Real-time detection of single electron tunneling through a T-shaped double quantum dot is sim- 
ulated, based on a Monte Carlo scheme. The double dot is embedded in a dissipative environment, 
and the presence of electrons on the double dot is detected with a nearby quantum point contact. 
We demonstrate directly the bunching behavior in electron transport, which leads eventually to a 
super-Poissonian noise. Particularly, in the context of full counting statistics, we investigate the 
essential difference between the dephasing mechanisms induced by the quantum point contact detec- 
tion and the coupling to the external phonon bath. A number of intriguing noise features associated 
with various transport mechanisms are revealed. 
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I. INTRODUCTION 

To control and manipulate electronic dynamics in 
nanoscale devices it requires knowledge of the involv- 
ing transport processes at single-electron level. The 
spectrum of current fluctuations, which characterizes the 
degree of correlation between charge transport events, 
serves as an essential tool superior to the average cur- 
rent in distinguishing various transport mechanisms. 1 ' 2 
Full counting statistics (FCS) of current^ has also been 
measured) 5 - - — owing in particular to the development of 
highly sensitive on-chip detection of single-electron tun- 
neling technique. All statistical cumulants of the number 
of transferred particles can now be extracted experimen- 
tally. 

Current fluctuations would obey a Poissonian process 
if the tunneling events were statistically independent. 
However, non-Poissonian fluctuation is in general a real- 
ity. In the case of transport through a localized state, the 
Pauli exclusion principle suppresses the noise^ leading 
to a sub-Poissonian statistics. Systems of multiple non- 
local states such as coupled quantum dots^ 2 - are more 
interesting. The intrinsic quantum coherence and many- 
particle interactions there result in different sources of 
correlations. The fascinating super-Poissionian noise 
thus occurs in various contexts, and thereby has been 
attracting a wide interest recently. 

A representing system, which will be studied in this 
work, is a T-shaped double quantum dot (TDQD) 
system^ - — as schematically shown in Fig.[TJ The sys- 
tem is of particular interest, as it can be mapped to a 
structure of quantum well in presence of an impurity in- 
side, which has been investigated experimentally*^ The 
source and drain electrodes of the TDQD are in such a 
configuration that maximizes locality versus nonlocality 
contrast. In addition, the TDQD is also influenced by an 
inevitable dissipative environment. The nearby quantum 




FIG. 1: Schematics of electron transport through a TDQD 
system, monitored continuously by the QPC current that de- 
pends on the charge state for excess electrons on two dots 
individually. The TDQD is also coupled to an inevitable dis- 
sipative phonon environment (not shown explicitly). 



point contact (QPC) that serves as a charge detector is 
asymmetrically coupled to the dots. The current through 
the QPC depends on the charge state, (rti,n 2 ), for excess 
electrons on two dots individually. Electron tunneling 
through the TDQD results in temporal changes in the 
charge state, and thus fluctuations in the QPC current. 

In this work, we first simulate the real-time detection 
of single electron tunneling through TDQD in experi- 
mental setup of Fig.[TJ The real-time simulation will be 
carried out based on a Monte Carlo approach to quantum 
master equation^ In particular, we demonstrate trans- 
parently the electron bunching, i.e., the super-Poissonian 
noise of charge transfer. The FCS is then studied, and the 
result shows good agreement with those obtained from 
Monte Carlo simulation. It therefore verifies the validity 
of our Monte Carlo method for real-time measurement. 
Particularly, it is found that the FCS is capable of dif- 
ferentiating between the dephasing effect induced by the 
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QPC charge detection and that by the coupling to the 
dissipative bath environment. A number of intriguing 
noise features will be revealed in association with vari- 
ous mechanisms, demonstrating thus the sensitivity and 
selectivity of cumulants to different transport properties. 

The paper is organized as follows. We describe the 
TDQD system and conditional quantum master equa- 
tion in Sec.|nJ The Monte Carlo simulation of real-time 
detection of single electron tunneling is summarized in 
Sec. lIII| together with the electron bunching behavior 
demonstrated vividly for the system in study. Having 
the FCS theory reviewed in Sec. lIVl we present the re- 
sults and discussions in Sec.|Vl Finally, we conclude in 
Sec.EH 



II. QUANTUM MASTER EQUATION THEORY 
OF QUANTUM MEASUREMENT 

A. Model device 

Let us start with the model setup for the QPC de- 
tection of single-electron tunneling through TDQD, as 
sketched in Fig.[T] The total Hamiltonian consists of the 
coupled dots system, the environment, and the coupling 
between them; i.e. H = H sys + H cnv + H sys ^ cnv . The 
electron Hamiltonian for the coupled dots system reads 

H sys = E (\ e Q^+^Q^cr)+'^2 Uim-fhii+U'hx^, (1) 
ct i 

where Q ZIJ = d\ a d lr7 - dl a d 2a , Qxa = d\ a d 2l j + d\ a d llT , 

ll i = So-^o-i ancl fycr = dialler, with di a {d\ a ) the elec- 
tron annihilation (creation) operator in the QD1 (I = 1) 
or QD2 (I = 2). Each quantum dot is assumed to have 
only one spin-degenerate level (a = f or \) in the bias 
window. The level detuning between the two dots is 
e = Ei — £2- The interdot coupling strength is f2. The 
intradot and interdot Coulomb interactions, Ui and U', 
are both assumed to be much larger than the Fermi lev- 
els. We shall be interested in the double-dot Coulomb 
blockade regime ; 18 ' 19 i.e. at most one electron resides in 
the TDQD. It can be realized in experiments by proper 
tuning the gate and bias voltages^ 

The environment is of H env = h p h + ^ a h a + /iqpc- 
It is composed of the phonon bath, the electron reser- 
voirs of the source and drain (a — L and R) elec- 
trodes, as well as the QPC detector. Each of them is 
modeled as a collection of noninteracting particles. The 
phonon bath assumes h p h — \ J2j f^j{p]+x 2 ). The elec- 
tron reservoirs are modeled with h a = ^ fe a t a kc} aktJ c a k a 

and h QPC = Y, P , a e v c U c P° + ^q,o- £q c \v c i°> respectively. 
These electronic parts are written in terms of electron 
creation and annihilation operators in the a-electrode 
and the QPC reservoirs states. 



The system-environment coupling can be written as 

a a, a 

+ E ^ QPC - ( 2 ) 

s=0 : l,2 

The first term describes the coupling with phonon bath, 
in which F p h = . XjXj. This term is responsible for the 
dot level energy fluctuations. The effect of phonon bath 
on the double-dot system is characterized by the phonon 
interaction spectral density, J p h(w) = l^il — ojj). 
The second term describes the transfer coupling between 
QD1 and the leads, in which f aa = J^k^akaCaka-- The 
third term describes the interaction between the QPC 
detector and the measured system, in which F^ PC = 
q aitspqcjxjCqa+H.c.). The amplitude of electron tun- 
neling through the QPC depends on the TDQD charge 
state h s — \s)(s\, with s = for no excess electrons, 
and s — 1 or 2 for one excess electron in QD1 or QD2, 
respectively. In other words, the current through the 
QPC sensitively depends on the charge states of the 
TDQD, and thus can be used to measure single elec- 
tron tunneling events. The effects of these electron reser- 
voirs components on the double-dot system are charac- 
terized individually by their interaction spectral densi- 
ties: J a (uj) = Y,ka \taka\ 2 S{uj - e ak ), and J s QPC (w) = 
E P ,,,a \t sp q\ 2 S(u - e pa + e qa ), respectively. 

In the Coulomb blockade regime, together with large 
source-drain voltage, the tunneling between QD1 and 
the electrode a =L or R can be characterized by the rate 
T a (ui) = 27r^ fco . \t a ka\ 2 &{w — e Q fc). In what follows we 
adopt flat bands in the electrodes, which yields energy in- 
dependent couplings Y a . Analogously, the coupling with 
the QPC is described by the rate T s = "2TTg p g q \t spq \ 2 VQPc, 
where Vqpc is the QPC bias voltage, g p and g q are the 
density of states in the QPC reservoirs. We assume here- 
after the density of states to be constant, and t spq reser- 
voir states independent. Thus T s just depends on the 
system charge occupation state |s)(s|, with s = being 
for zero excess electron, and s = 1 or 2 for one excess 
electron in QD1 or QD2, respectively. For the QPC cur- 
rent, I s = eT s , we have Iq > 1% > I±, as implied in the 
scheme of Fig.[TJ 



B. Quantum master equation theory 

To describe the quantum measurement, we exploit 
the reduced density operator p^ n '(t) of the TDQD sys- 
tem, for the specified number n electrons having passed 
through the QPC detector and being recorded up to 
the given time. The related conditional quantum mas- 
ter equation can be derived, which is greatly simplified 
under the Born-Markov approximation ! 18 ' 19 ' 21 ' 22 Here, 
instead of using the "n" -resolved equation directly, it 
is convenient to introduce its x _s P ace counterpart via 
p{x,t) = ^2 n e mx p^ n \t), where x is tne so-called count- 
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ing field. Let the quantum master equation be formally 
^- t p( X ,t) = -(hC + 1l x )p(x,t), (3) 

where C ■ = [H sys , ■ ] is the system Liouvillian, and 1Z X is 
the dissipation superoperator to be specified soon. Here- 
after, we set unit of h = e = 1 for the Planck constant 
and electron charge, unless where is needed for clarity. 

To that end, let us recast the reduced density matrix 
in the vector notation, 

P = (POO, Pll, P22, Pl2, P21, Pll, P22, Pl2, P2lf ', (4) 

with p ss , = (sf|p|st) and p ss > = (sl\p\sl). There 



are 9 nonzero elements of the reduced density matrix in 
the Coulomb blockade regime. The other 16 elements, 
describing coherence between different spin states, are all 
zeroes as the system-environment coupling consider here 
does not cause spin flip. The dissipation superoperator 
1Z X in Eq. ([3]) is simply a 9 x 9 matrix. 

In the large voltage limit and Coulomb blockade 
regime, the involving Fermi functions can be approxi- 
mated by the step function of either one or zero. Un- 
der the Born-Markov approximation for weak tunnel 
coupling and second-order perturbation theory in the 
electron-phonon coupling, the dissipation superoperator 
1Z X matrix is obtained explicitly 



/2r L + r g(x) 


-r R 











-r R 








\ 


-r L r R 































T 2 q(x) 























r+ 


r_ 


r d 




















r+ 


r_ 





r d 














-r L 











o r R 


f7I«(x) 





























T2q{x) 























r+ 


r_ 


r d 





V o 














r+ 


r_ 








Here we have introduced 



1± 



4-coth f A ) 
A \2knTJ 



and 



T d = 2( r R + 7 P h + 7 X ): 



il 2 /A 
7 P h = 2tt— J ph (A)coth f — T . , 



7x = (V7i 



2 v / 7T7f<7(x), 



(6) 



(7) 



(8) 



(9) 



with q(x) = 1 — e' x , A = ^ e 2 + 4S1 2 the Rabi frequency 
of the double-dot system, fc R the Boltzman constant, and 
T the temperature. The physical processes described by 
the transfer elements in Eq. §5§ are clearer— 

Consider 1Z\ lm = — T L , 7£q n = — r R , and their op- 
posite spin counterparts, as inferred from Eq. ([5]). These 
identities agree with the sequential tunneling picture, 
with Tl being the rate of electron tunneling from the 
left electrode to system, and T R being that from system 
to the right electrode; both via QD1. 

The parameter T± [Eq. ©] denotes the nonsecular el- 
ements for the population-to-coherence transfers: T + = 



^12,11 = ^21,11 ) r - = ^12,22 = ^21,22: and their °PP°- 

site spin counterparts, as denoted in Eq. ([S]). They are 
purely due to the phonon bath environment. 

The parameter [Eq. ([7])] is the total decoherence 
rate between two levels: Td = TZ\ 2 12 = ^21 21 and ^he 
opposite spin counterparts, as denoted in Eq. (J3J) - The 
total decoherence rate is composed of not just T R due 
to the electron depopulation to collector, but also 7 P h 
[Eq. flSJ] and j x [Eq. ©], due to the phonon bath cou- 
pling and the QPC detection, respectively. Note that 
7d = 7 x =o = (V7T — n/TI) 2 denotes the dephasing rate 
induced by the existence of QPC ensemble. 

In the numerical demonstrations below, we set the 
phonon bath spectral density Jph(w) = 2?7we~ w / Wc . Here, 
the dimensionless parameter 77 reflects the strength of dis- 
sipation and uj c is the Ohmic high energy cutoff. 



III. MONTE CARLO SIMULATION OF SINGLE 
ELECTRON TUNNELING 

Experimentally, the most intuitive method for measur- 
ing the FCS of electron transport is to count electrons 
passing one by one through the conductor. The real- 
time detection of single electron transport enables direct 
evaluation of the probability distribution function of the 
number of electrons transferred through device within a 
given time period. In addition to the current and the 
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FIG. 2: Time traces of the QPC current fluctuations correspond to different charge states in the TDQD (shown on the right). 
The arrows indicate transitions where an electron is entering the QD1 from the left lead, (a) and (b) correspond to different 
temperatures T = fOOmK and T = 800 mK, respectively. Other parameters are: Tr = 4l\ = Q = 20 kHz, e = —40, 
t) = 2 x 10 -6 , and ut c = 0.4 meV. The time step used is r=0. 002ms, such that the minimum characteristic time-scale of the 
system can be clearly resolved. 



shot noise, which are the first and second moments of 
this distribution, this method gives also access to higher 
order cumulants. 

Here we outline the Monte Carlo simulation scheme 
for single electron tunneling through deviceJ^ Let us 
start with the quantum master equation ([3]) and denote 
C x = — iC — 1Z X there. The solution to the reduced state 
simply reads p(x, t) — e CxSt p(x, to), for an arbitrary ini- 
tial condition p(x,to) = p{to) and finite time interval 
St = t — to- Ts counterpart in the particle number "n"- 
space is obtained via the inverse Fourier transform 

2 ^X 



,C x St-inx 



2tt 



p(t ) = U(n,5t)p(t ). (10) 



The involving propagator U(n,dt) is completely deter- 
mined by the dynamic structure of the master equation 
([3]), regardless of the initial state. Therefore, we can nu- 
merically evaluate it by a "one-time task", such as fast 
Fourier transform, which results in an efficient real-time 
simulation. 



electrons having passed through QPC during [tj,tj +r]. 
The second equality of Eq. (fTTj) implies that if we stochas- 
tically generate rij according to Pr (rij) for each time in- 
terval [tj , tj + t] , and collapse the state definitely onto 
p c (rij,tj + t), we have in fact simulated a particular re- 
alization for the selective state evolution conditioned on 
the specific measurement results. 

For the output current in a particular real-time mea- 
surement, we have 

W*) =hp00 + h(pll+ph) + l2(p C 22 + Pk)+Z(t). (12) 

The first three terms determine the conditional evolution 
of the charge state. The last term £ (t) originates from the 
intrinsic noise of detector. Here, we consider in the diffu- 
sive regime, where £(t) is a Gaussian variable with zero 
mean value and the spectral density S% = 2e(/ QPC ), with 
(/ QPC ) the average stationary QPC current. Accordingly, 
we can stochastically generate rij, the number of elec- 
trons having passed through QPC during [tj,tj + r], via 



Specifically, consider the evolution of state p(tj) at tj n i ~ /*/ dt>I QPc(t') - [IoPoo + h(pu + Pii) + h{ph + 
to p^itj+r) at tj+r: p^(tj+r) = U(rij,T)p(tj), and Ph)) T + where dW(tj) is the Wiener increment 



denote Pr(n 3 ) = Tr[p^ n ^(tj + r)] that is the probability 
of having rij electrons passed through QPC during the 



during [tj , tj + r] . 

A typical example of simulated real-time detector cur- 



time interval [tj , tj + r] . If the measurement is made but rent is displayed in Fig.rj for (a) T - 100 mK and 



the result is ignored, then the (mixture) state reads 



p(tj 



(ii) 

Here, p c (nj,tj + r) = p^\tj + r)/Pr(nj) is the nor- 
malized state, conditioned by the definite number of rij 



(b) T — 800 mK, respectively. The temperatures are 
much smaller than the bias voltage across the TDQD, as 
well as the interdot and the intradot charging energies. 
Thus, electrons are transported in one direction, i.e. from 
the left lead to the right one. The corresponding charge 
states (ni,ri2) of the TDQD, with ri\ and ni being the 
excess electrons in QD1 and QD2, respectively, are shown 
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FIG. 3: Statistical distribution of the number N of electrons 
entering the TDQD during a given measurement time span 
of t c = 3 ms for (a) T = 100 mK, and (b) T = 800 mK, 
respectively. The other parameters are the same as in Fig.[2] 
The solid line shows the distribution calculated from Eq. {16} . 



on the right. We choose rates of tunneling to the left and 
right leads as IY = 5 kHz and Tr = 20 kHz, such that all 
the tunneling rates are within the bandwidth (~ 30 kHz) 
of the QPC detector^ The QPC conductance (without 
excess electrons in TDQD) is G = 0.02 e 2 /h, and QPC 
bias voltage V'qpc = 0.15 mV, which corresponds to a 
QPC current Jo « 0.12 nA. The conductance is assumed 
to decrease by 4% or 2%, when an electron occupies on 
the QD1 or QD2, respectively. The corresponding QPC 
currents are I\ ~ 0.115 nA and I2 ~ 0.118 nA. The rate 
of dephasing due to QPC charge detection- 6 is then de- 
termined as 7d = -(y/Ti — VI2) 2 ~ 75kHz. 

In Fig. [21 each arrow indicates the tunneling of one elec- 
tron from the left lead into QD1. Accordingly, a charge 
state transition (0,0)— >(1,0) occurs. Remarkably, we ob- 
serve a clear signature of electron bunching, especially 
when temperature is high, as shown in Fig.^b). This 
phenomenon can be understood in terms of dynamical 
channel blockade, as will be explained later in more de- 
tail. For a thorough analysis of the statistical properties 
for electron transport through the TDQD device, we re- 
sort to the FCS. It is enabled by the probability distri- 
bution for the number N of electrons entering the QD1 
from the left lead (the down steps indicated by arrows 
in Fig. [2]) during a given time span £ c , by counting the 
number of electrons from the time traces analogous to 
the ones shown in Fig. [5] The obtained current distribu- 
tions for T — 100 mK and T — 800 mK are displayed by 
histograms in Fig.|3fa) and (b), respectively, by repeat- 
ing this counting procedure on one thousand independent 
traces with the measurement time span of t c = 3 ms. 



IV. FULL COUNTING STATISTICS THEORY 

We consider hereafter the counting statistics of cur- 
rent through the TDQD device. The central quantity 
is then the probability that a given N electrons have 
passed through the device during the counting measure- 
ment time span t c . This probability is related to the 
particle-number-resolved reduced density operator as 



P(N, t c ) = Tip( N \t c ), where the trace is over the system 
degrees of freedom. Note here "JV" denotes the number 
of electrons transferred through the device, rather than 
the QPC detector considered earlier. From the knowl- 
edge of these probabilities one can easily derive not only 
the current and noise, but all the cumulants of the cur- 
rent distribution. The associated cumulant generating 
function (CGF) g(ip) reads 

efffoO ^J2 p ( N ^c)e~ iNv , (13) 

N 

where ip is the counting held on a specified TDQD lead. 
All cumulants of the current can be obtained from the 
CGF by performing derivatives with respect to the count- 
ing held 

(I k ) = -{-i dtp ) k g{y)\ v = Q . (14) 

The hrst three cumulants are related to the average cur- 
rent, the (zero-frequency) current noise, and the skew- 
ness, respectively. 

To evaluate the CGF, let us consider g((p, t) = 
P^ N \t)e lNv . Its equation of motion reads 

Q{(p,t)=C v g{(p,t). (15) 

The involving generator C v is completely determined 
by the associated conditional master equation,!? which 
is similar to Eq. ([3]) but with the number of electrons 
passing through the TDQD device being resolved. The 
formal solution to Eq. (fT5j) is g((p,t) = e Clfit g(ip, 0). 
Straightforwardly, the CGF is determined as g(<p) = 
— ln{Trg(ip, t c )}. Actually, we are most interested in the 
zero- frequency limit, i.e. the counting time t c is much 
longer than the time of tunneling through the system. 
The CGF then simplifies to g(<p) = — A m ; n (<^)t c , where 
A m i n (</ ? ) is the minimal eigenvalue of C v that satisfies 
A mi n(y ->())-► 0i^2a With the knowledge of CGF, 
the distribution function can be readily obtained via£ 

P(jV) = *e. e -n(.<P)-Uf<P, (i 6 ) 
Jo 27r 

V. RESULTS AND DISCUSSIONS 

In what follows, we focus our analysis on electron tun- 
neling from the left lead to QD1. The relevant N is 
then the number of electrons entering the dot from left 
lead, and <p is the corresponding counting held. Note the 
same calculations apply to the counting for the right lead. 
The numerical results for the probability distribution are 
plotted by the solid lines in Fig. [31 It shows a striking 
agreement with the histograms obtained by Monte Carlo 
simulation; thus, it also verifies the validity of our Monte 
Carlo method for real-time measurement. The two dis- 
tributions in Fig.[3{a) and (b) are rather different. The 
latter shows a broader and more asymmetric distribution 
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than the former. We characterize the differences quanti- 
tatively based on the FCS analysis, as follows. 

Consider first the situation without electron-phonon 
interaction (rj = 0). The first current cumulant gives the 
average current (I) = 2r L r R /r cff , with r cff = 4r L + T R 
the total effective tunneling widthi 18 i 19 It is independent 
of level detuning e, interdot coupling ft, and QPC charge 
detection induced dephasing 7^. In contrast to the cur- 
rent, valuable information can be extracted in the second 
cumulant (zero-frequency shot noise). By expressing it in 
terms of the Fano factor F = (I 2 ) / (I), we readily obtain 

p , sr L r R 2 (7d + r R )2 + 4e 2 

b ~ l ^2 r 21 l! R -t . 2 „ ■ [If) 

r cs (7d + r R )r^ ff ft 2 

It thus allows us to get more knowledge than the current 
on the processes involved in the electronic transport. In 
Fig.0Ja), the Fano factor F is plotted against the level 
detuning e for different QPC-induced dephasing rates 7d- 

Noticeably, a significantly enhanced Fano factor is ex- 
pected when the dot levels are far from resonance. In this 
case, electron transitions between the two dots are sup- 
pressed. For instance, if an electron is tunneled into QD2, 
it will dwell on it for a long time. In the strong Coulomb 
blockade regime, the electron in QD2 will block the cur- 
rent until it is removed and tunneled out to the right 
lead. Consequently, a mechanism of dynamical channel 
blockade is developed ) 19 ' 24 i 33 ~ — and the system exhibits 
strong electron bunching behavior, which leads eventu- 
ally to a profound super-Poissonian noise. 

It is worth noting that a strong super-Poissonian noise 
is more readily achieved for a small Jl [cf. Eq. (fT7|) ]. In 
this case, electron tunneling between the two dots is 
suppressed, which would enhance electron localization of 
electron in QD2. Electrons thus tend to be transferred 
in bunches. Remarkably, as SI — > a diverging Fano fac- 
tor is found. Here, the strong bunching behavior is re- 
sponsible for the divergency^ Similar results were also 
reported in Ref. |H, where counting statistics for electron 
transport through a double-dot Aharonov-Bohm inter- 
ferometer was investigated. However, in that case the di- 
vergence is closely related to a separation of the Hilbert 
space of the double-dot into disconnected subspaces that 
contain the spin singlet and triplet states for double oc- 
cupancy. 

The QPC charge detection has a two-fold effect on the 
measured system. On one hand, it causes dephasing be- 
tween the two dot-levels. Generally, the dephasing mech- 
anism leads to the suppression of noise^ which explains 
the 7d-dependence of the Fano factor shown in Fig.0] 
i.e. the noise is reduced with rising dephasing rate, par- 
ticularly in the regime far from resonance. On the other 
hand, the measurement gives rise to the so-called quan- 
tum "Zeno" effectj22r— which dominates in the regime 
of large dephasing rate, and results in a strong dynamic 
charge blockade behavior. The noise is finally enhanced 
with increasing dephasing rate, as we have checked (not 
shown in Fig.@|. 
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FIG. 4: (a) Fano factor and (b) normalized skewness ver- 
sus level detuning for different QPC-induced dephasing rates. 
The electron charge is set to be e = 1. The temperature is 
T=100mK, and the other parameters are the same as in Fig. [2] 
but in the absence of phonon bath; i.e., rj = 0. 



The Fano factor studied so far proves to be much sen- 
sitive than the average current, it only reveals, however, 
limited information about the QPC charge detection in- 
duced dephasing mechanism. We thus expect more infor- 
mation to be extracted in the next order cumulant, i.e. 
the skewness. The numerical results for the normalized 
skewness S = (I 3 ) / (I) is displayed in Fig.|U(b) as a func- 
tion of level detuning e. Without QPC charge detection 
induced dephasing, the skewness reaches the maximum 
when the dot levels are in resonance (e = 0). As the de- 
phasing rate grows, it turns into a local minimum. At the 
edges of the resonance a double maximum structure sym- 
metric around e = is observed. The local maxima shift 
away from resonance with increasing dephasing rates, as 
clearly demonstrated in Fig.HJb). 

Now let us turn to the influence of phonon heat bath 
that induces dephasing between two dots. We will reveal 
the essential difference between the dephasing induced 
via QPC charge detection and that by phonon coupling. 
The former can be modified via the coupling between the 
TDQD and the QPC, while the latter is generated with 
emission and absorption of phonons and increases with 
rising temperature. Here, we limit our discussions to the 
temperatures well below the Coulomb charging energies 
and the bias voltage. Therefore the temperature acts 
solely due to the coupling to the phonon heat bath. 

The calculated Fano factor and normalized skewness 
versus level detuning are plotted in Fig. [5] for different 
temperatures. Increasing temperature results in the en- 
hancement of phonon-bath induced dephasing, as the 
emission and absorption of phonons occur more fre- 
quently. Analogous to the QPC detection, it gives rise 
to a mechanism of localization in QD2, and consequently 
to a dynamical channel blockade. Eventually, electron 
transport through QD1 occurs in bunches during lapses 
of time when the QD2 is empty. This is confirmed by 
the real-time trajectory [cf. Fig.^b)]. Thus the noise in- 
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FIG. 5: (a) Fano factor and (b) normalized skewness versus 
level detuning for different temperatures with 7d = 75 kHz. 
The electron charge is set to be e = 1. The other parameters 
are the same as in Fig. [2] including those of phonon heat bath 
of rj = 2 x 10" 6 and lj c = 0.4 meV. 



creases with rising temperature, as displayed in Fig.JSJa). 

If the TDQD is coupled only to the QPC, the cumu- 
lants are symmetric around e = 0, as shown in Fig. 2] 
However, with non-zero coupling to the heat bath the 
noise exhibits a clear asymmetry (see Fig. [5]). This is 
essentially due to another consequences of phonon cou- 
pling: The phonon mediated transition can partially re- 
solve the dynamical charge blockadei 24 i 40 For instance, at 
zero temperature and e > 0, spontaneous phonon emis- 
sion can partially lift the localization caused by dephas- 
ing, and the spectrum turns out to be asymmetric. Simi- 
lar argument applies to the regime of finite temperature, 
where both phonon absorption and emission take place. 
The fact that emission is more likely than absorption ex- 
plains the observed asymmetry 

The skewness of the current distribution, which was 
not explored in Ref. 24, is found to be more sensitive to 
thermal phonon bath-induced dephasing. Without cou- 
pling to the phonon bath, the skewness shows a symmet- 
ric double maximum structure. The spectral becomes 
asymmetric in the presence of electron-phonon interac- 
tion. With increasing temperature, the phonon absorp- 
tion takes place more frequently at e < , and thus 



greatly enhance the local maximum. While in the oppo- 
site regime of e > the local maximum is reduced due to 
suppressed phonon emission. 



VI. SUMMARY 

We have investigated electron transport through a T- 
shaped double quantum dot system by utilizing a quan- 
tum master equation approach. The major advantage of 
the present approach is its simplicity of treating prop- 
erly the dephasing mechanism of the QPC charge detec- 
tion and that due to an external heat bath. In addition, 
this approach has the merit of dealing with other sources 
of dephasing, such as that entailed by anti-resonances. 41 
Particularly, based on the Monte Carlo scheme, real-time 
detection of single electron tunneling is simulated by ex- 
ploiting the sensitivity of a current, passing through a 
nearby quantum point contact, to the fluctuating charge 
on the quantum dots. Owing to the interplay between 
the Coulomb interactions and the dephasing mechanisms, 
a strong bunching behavior in the charge transfer was 
detected, which leads eventually to a super-Poissonian 
noise. 

Furthermore, full counting statistics of the transport 
current is analyzed based on the probability distribu- 
tion, which is determined by ensemble average over a 
large number of single trajectories. It is demonstrated 
that the dephasing mechanism of the QPC charge detec- 
tion and that owing to the external heat bath give rise 
to distinct and intriguing features. It thereby enables 
us to achieve a clear identification of different dephas- 
ing sources. Investigations of various processes involved 
in the electronic transport through similar devices are 
highly desirable in experiments. 



Acknowledgments 

Support from the National Natural Science Foundation 
of China under Grants No. 10904128, and the Research 
Grants Council of the Hong Kong Government (Grant 
No. 604709) are gratefully acknowledged. 



* Electronic address: jyluo@zust.edu.cn 
t Electronic address: yyan@ust.hk 

1 Y. M. Blanter and M. Biittiker, Phys. Rep. 336, 1 (2000). 

2 Y. V. Nazarov, Quantum Noise in Mesoscopic Physics 
(Kluwer, Dordrecht, 2003). 

3 L. S. Levitov, H. W. Lee, and G. B. Lesovik, J. Math. 
Phys. 37, 4845 (1996). 

4 D. A. Bagrets and Y. V. Nazarov, Phys. Rev. B 67, 085316 
(2003). 

5 W. Lu, Z. Ji, L. Pfeiffer, K. W. West, and A. J. Rimberg, 
Nature 423, 422 (2003). 

6 J. Bylander, T. Duty, and P. Delsing, Nature 434, 361 



(2005). 

7 R. Schleser, E. Ruh, T. Ihn, K. Ensslin, D. C. Driscoll, and 
A. C. Gossard, Appl. Phys. Lett. 85, 2005 (2004). 

8 L. M. K. Vandersypen, J. M. Elzerman, R. N. Schouten, 
L. H. W. van Beveren, R. Hanson, and L. P. Kouwenhoven, 
Appl. Phys. Lett. 85, 4394 (2004). 

9 S. Gustavsson, R. Leturcq, B. Simovic, R. Schleser, T. Ihn, 
P. Studerus, K. Ensslin, D. C. Driscoll, and A. C. Gossard, 
Phys. Rev. Lett. 96, 076605 (2006). 

T. Fujisawa, T. Hayashi, R. Tomita, and Y. Hirayama, 
Science 312, 1634 (2006). 

1 L. Y. Chen and C. S. Ting, Phys. Rev. B 43, 4534 (1991). 



8 



13 



12 W. G. van der Wiel, S. D. Franceschi, J. M. Elzerman, 
T. Fujisawa, S. Tarucha, and L. P. Kouwenhoven, Rev. 
Mod. Phys. 75, 1 (2003). 

T.-S. Kimand and S. Hershfield, Phys. Rev. B 63, 245326 
(2001). 

14 P. S. Cornaglia and D. R. Grempel, Phys. Rev. B 71, 
075305 (2005). 

15 I. Djuric, B. Dong, and H. L. Cui, Appl. Phys. Lett. 87, 
032105 (2005). 

16 A. Nauen, I. Hapke-Wurst, F. Hohls, U. Zeitler, R. J. Haug, 
and K. Pierz, Phys. Rev. B 66, 161303 (2002). 

17 S.-K. Wang, J. S. Jin, and X.-Q. Li, Phys. Rev. B 75, 
155304 (2007). 

18 J. Y. Luo, X.-Q. Li, and Y. J. Yan, Phys. Rev. B 76, 
085325 (2007). 

19 J. Y. Luo, X.-Q. Li, and Y. J. Yan, J. Phys.: Cond. Matt. 
20, 345215 (2008). 

20 K. Ono, D. G. Austing, Y. Tokura, and S. Tarucha, Science 
297, 1313 (2002). 

21 X. Q. Li, P. Cui, and Y. J. Yan, Phys. Rev. Lett. 94, 
066803 (2005). 

22 X. Q. Li, J. Y. Luo, Y. G. Yang, P. Cui, and Y. J. Yan, 
Phys. Rev. B 71, 205304 (2005). 

23 T. Brandes, Phys. Rep. 408, 315 (2005). 

24 G. Kiefilich, E. Scholl, T. Brandes, F. Hohls, and R. J. 
Haug, Phys. Rev. Lett. 99, 206602 (2007). 

25 R. Aguado and T. Brandes, Phys. Rev. Lett. 92, 206601 
(2004). 

26 S. A. Gurvitz, Phys. Rev. B 56, 15215 (1997). 



27 C. W. Groth, B. Michaelis, and C. W. J. Beenakker, Phys. 
Rev. B 74, 125315 (2006). 

28 C. Flindt, T. Novotny, and A.-P. Jauhol, Europhys. Lett. 
69, 475 (2005). 

29 G. Kiefilich, P. Samuelsson, A. Wacker, and E. Scholl, 
Phys. Rev. B 73, 033312 (2006). 

30 A. N. Korotkov and D. V. Averin, Phys. Rev. B 64, 165310 
(2001). 

31 S. A. Gurvitz, L. Fedichkin, D. Mozyrsky, and G. P. 
Berman, Phys. Rev. Lett. 91, 066801 (2003). 

32 J. Y. Luo, H. J. Jiao, F. Li, X.-Q. Li, and Y. J. Yan, J. 
Phys.: Cond. Matt. 21, 385801 (2009). 

33 A. Cottet, W. Belzig, and C. Bruder, Phys. Rev. Lett. 92, 
206801 (2004). 

34 M. Gattobigio, G. Iannaccone, and M. Macucci, Phys. Rev. 
B 65, 115337 (2002). 

35 R. Sanchez, G. Platero, and T. Brandes, Phys. Rev. Lett. 
98, 146805 (2007). 

36 G. Kiefilich, H. Sprekeler, and E. Scholl, Semicond. Sci. 
Technol. 19, S37 (2004). 

37 A. Cottet and W. Belzig, Europhys. Lett. 66, 405 (2004). 

38 F. Li, H. J. Jiao, J. Y. Luo, X.-Q. Li, and S. A. Gurvitz, 
Physica E 41, 1707 (2009). 

39 D. Urban and J. Konig, Phys. Rev. B 79, 165319 (2009). 

40 A. Braggio, C. Flindt, and T. Novotny, J. Stat. Mech.: 
Theory Exp. 2009, P01048. 

41 L. E. F. F. Torres, H. M. Pastawski, and E. Medina, Eu- 
rophys. Lett. 73, 164 (2006). 



